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Abstract 

Magnetic susceptibility measured for CaA^Og is analyzed by a 
method with high temperature expansion. The analysis is consistent 
with the 5-value which is close to 2 and is observed by an ESR ex- 
periment. Four exchange parameters in a two-dimensional Heisenberg 
model are determined to represent CaV^Og. All the exchange param- 
eters are about 500 K and strongly frustrated. The observed spin gap 
originates from the frustration. 

keywords: CaV 4 Og, spin gap, magnetic susceptibility, high temperature ex- 
pansion, two-dimensional antiferromagnet, frustration 



1 Introduction 

Spin gaps in low dimensional materials have been investigated by a num- 
ber of researchers. While origins and effects of spin gaps are fundamentally 
and generally interesting, those for two-dimensional materials are especially 
important in possible relevance to the high temperature superconductivity. 
Recently Taniguchi et al. U found a spin gap in CaV40g by magnetic suscep- 
tibility and NMR measurments. CaV40g is a layered insulator and has spin 
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Figure 1: Lattice structure for vanadium spins in a layer of CaX^Og. The 
Heisenberg model examined in text includes 4 dependent exchange parame- 
ters J e (bold solid line), J' e (solid line), J c (bold dashed line) and J' c (dashed 
line) . 



degrees of freedom at vanadium ions. || This is the first clear observasion of 
a spin gap in a two-dimensional spin system. 

Theoretical efforts have been done to describe CaA^Og as a two-dimensional 
Heisenberg model and to explain the origin of the spin gap. H f|, ||, 0, |[ 
H, [10], ?, [12|, [13], |14| Almost models considered are included, as special cases, 



in a Heisenberg model with 4 kinds of exchange parameters; J e for edge- 
sharing plaquette links, J' e for edge-sharing dimer links, J c for corner-sharing 
plaquette links and J' c for corner-sharing dimer links as shown in Fig. [l]. In 
these studies authors assumed values or constraints for the parameters and 
constructed their theories. Among them, Troyer || performed a quantum 
Monte Carlo simulation and argued that the system does not form a spin gap 
in the unfrustrated homogenous case of J e = J' e ^ and J c = J' c = 0. The 
spin gap is possibly formed by inhomogeneity: i. e. the difference between 
J e and J' e . It is also possible that the spin gap originates from frustration 
due to nonzero J c and J' c . 

In a previous paper, we roughly estimated J e and J c on the condition 
of J' e = J e and J' c = J c . We fitted a high-temperature expansion to the order 
of 1/T 3 for the susceptibility to experimental data. The result was J e ~ 
610 K and J c ~ 150 K with keeping the opposite possibility. |L5| Calculation 
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of numerical diagonalization shows that the frustration among interactions 
with J e and J c enhances the spin gap in comparison with the case of J c = 0. 
This enhancement agrees with the result of a perturbation calculation. [|j] 
Most recent experiments give more detailed information on this material: 



i) Neutron scattering measurement |16) suggests J c >J e . (ii) ESR measure 



ment |T7[ precisely determined the g- value as g = 1.96, which is rather close 
to 2. The Curie constant is then C = 0.003713 emu/g. As for our previous 
paper, experiment (i) rather supports the opposite possibility for the values 
of the exchange parameters; i. e. J c ~ 610 K and J e ~ 150 K. On the 
other hand, we have used g = 2 of a free spin for g-value without explicit 
notice there. This g-value is close to the truth. Gelfand et al. estimated the 
exchange parameters as J e ~ J' e ~ 190 K and J c ~ J' c ~ J e /2 by various 
expansions. In the estimation, they determined the g- value as g = 1.77. The 
difference between this value and the experimental value seems to be hardly 
explained. 

In view of the above, we expect that the previous estimation for J c and 
J e is roughly realistic, if we exchange the values of them. However, the esti- 
mation still includes unsatisfactory points. First we examined a Heisenberg 
model including only two independent exchange parameters, J e and J c . It is 
at least necessary to independently treat parameters for plaquette and dimer 
links in order to determine the origin of the spin gap. Second we have a room 
to improve the fitting method. We obtained an expansion of the susceptibil- 
ity of the Heisenberg model up to a finite order and fitted it directly to the 
experimental data. This method is rather usual for high temperature fitting. 
However the estimation is rough unless fairly high order expansion is carried 
out, and it is hard to calculate high order terms of expansion for a system 
with many independent exchange parameters. 

In this paper, we determine the exchange parameters of the Heisenberg 
model to describe CaA^Og in a novel fitting method. Following this method, 
we first construct an experimental formula for the magnetic susceptibility in 
power series of 1/T. This procedure is explained in detail in section 2. On the 
other hand, we obtain a high temperature expansion of the susceptibility for 
a Heisenberg model with 4 independent exchange parameters; the expansion 
coefficients are functions of the exchange parameters. We determine the 
parameters by fitting the expansion coefficients to the coefficients of the 
experimental formula. The fitting is carried out in section 3. The result 
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shows that the origin of the spin gap for this material is the frustration 
rather than the inhomogeneity. Discussion is devoted in section 4. 



2 Experimental formula for susceptibility at 
high temperatures 

The experimentally obtained magnetic susceptibility, x EXP (T), is a set of 
data for T<700 K. We construct an experimental formula for the data in a 

power series of 1/T. The formula will be used to compare the experiment to a 
theory at high tempratures. To make expansion coefficients dimensionless we 
use x = Tq/T as the expansion parameter, where To is an arbitrary constant 
with the dimension of temperature. We hereafter employ To = 700 K without 
spoiling generality. The formula is then written as 

xCD = hm x {n \n 

fi n 

x (n) (T) = -[i + £A^ m ] (l) 

1 771=1 

with Curie constant C deteremined by the ESR measurement. Here the 



expansion coefficients v4 m 's are fitting parameters. JT| 

To determine An's we introduce a novel fitting method as will be ex- 
plained below. We first define the quantity 

oo 

<f> m {x) = A m + J2 A m+i x l (2) 
i=i 

for each m. It is important that the quantity reduces to coefficient A m in 
the high temperature limit: 

m (O) = A m . (3) 
4> m (x) is also defined by the recursion equation: 

T 

M x ) = 75 x(T), 

0mO) = (0 m _l(x) - An-l)- (4) 
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for m = 0,1, 2, • • • with Aq = 1. The experimental data corresponding to 
4> m {x) is similarly defined in the recursive transformation: 



Here m fits 4> m because \ is constructed to fit x ■ It is noticed that 
4> m is only formally introduced since the starting function \ is now unknown 
and is the object which we will obtain finally On the other hand we have 
the data x EXP an d so 4> m xp i s actually obtained as is mentioned below. 

Now we determine {A m } together with {4> m XP {x)} one by one. For m = 
0, we have (pQ XP (x) by multiplying T/C to x EXP (T) and have A = 1 by 
the definition. Next we obtain <pf xp (x) by the transformation (pf xp (x) = 
((pQ XP (x) — Aq)/x in eq. (fjD with m = 1. To obtain Ai we use a simple 
fitting function fi(x) = a\ + b\ exp(— c\x) and determine the parameters a\, 
b\ and c\ to make fi(x) fit 4>f xp (x) by the least square method. Using these 
parameters, A\ is given by A\ = /i(0) = a\ + b\ corresponding to A\ = <pi(0) 
in eq. (H) with m — 1. By repeating this process, we obtain (j> m {x) and v4 m for 
arbitrary m: When m _i(x) and v4 m _i have been known, <p^ l xp (x) is given 
by eq. (H). Then we make function 



fit 0^ xp (x) and determine parameters a m , b m and c m . A m is obtained by 
A m = /m(0) = a m + b m corresponding to eq. (||). Thus we can determine any 
coefficient A m inductively. 

In Fig. ||, we show 4> m xp (x) along with f m (x) for m — 1 to 4. Here we 
have used a weight function exp(— 1000/T) for fitting by the least square 
method. Optimal values for a m , b m , c m and A m are shown in Table |]. In 
principle we can obtain coefficient A m for any m. However, experimental 
data 4> m xp (x) becomes dispersive as m increases, so that A m with very large 
m cannot be obtained. When we change the weight function, the values of 
Am's change with a strong correlation: |A m |'s have a tendency to increase 
or decrease simultaneously. The error depending on the choice of the weight 
function seems to be roughly several percents. 

We have obtained values for A m 's within some accuracy. They are inde- 
pendent of n if n is sufficiently large. It is now instructive to substitute the 




^{X) = (^f(x)-An-l)". 



(5) 



fm\£) Q"m ^m6Xp( C m X) 



(6) 
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Figure 2: Transformed experimental data m xp (a;) and fitting function f m (x) 
for m — 1 to 4. 



T able 1: Fitting parameters and coefficients A m 's. 



m 


o m 




Cm 


A 


1 


-0.20907 


-0.76379 


1.2496 


-0.97286 


2 


0.15067 


0.77787 


0.68028 


0.92854 


3 


-0.10418 


-0.42303 


0.41366 


-0.52721 


4 


0.04797 


0.12828 


0.29549 


0.17624 


5 


-0.01365 


-0.02554 


0.26649 


-0.03918 


6 


0.00249 


0.00484 


0.26474 


0.00733 


7 


-0.00024 


-0.00104 


0.15695 


-0.00127 


8 


0.00002 


0.00014 


0.08627 


0.00016 
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Figure 3: X s as functions of T. They are given by eq. ([[]) with A m 's 
in Table [I]. Data of the experimental susceptibility x EXP are a l so shown by 
open circles. 

values for A^s into x i n ec L- (H) with small n as well as large n. Results 
for several values of n are shown in Fig. [3]. x approaches to x EXP as n 
increases, confirming the validity of this method. 

3 High temperature expansion for a Heisen- 
berg model and determination of exchange 
parameters 

We assume that spins at vanadium sites in a layer of CaX^Og are described 
by a two-dimensional Heisenberg model which is represented in Fig. |I[ The 
Hamiltonian is written as 

H = £ JyS, • Sj, (7) 

<i,j> 

where Si is the spin at site i. The exchange parameter is J e , J' e , J c or J' c 
if it corresponds to a link indicated in Fig. [l| and is zero otherwise. 

By the high temperature expansion, magnetic susceptibility of a Heisen- 
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berg model is generally written in the following form: 

X HTE {T) = '-[l+Y.F^ m ] (8) 

1 m=l 

with x = Tq/T and the experimentally determined Curie constant C. We 
have used T = 700 K as in the previous section. The coefficients F m 's are 
functions of the exchange parameters and are calculated by the standard 
diagramatic method. [IS 



In the present model of eq.(|7|), the coefficients F m 's are functions of J e , 
J' e , J c and J' c . We obtained them for m = 1, 2 and 3 as 

F X T Q = -l(2J e + J' e + 2J c +J' c ), 

F 2 T 2 = L[(AJ e - J'JJ' e + 2(2 J e + f e )(2J c + f c ) + (AJ c -f c )J' c }, 

= ^[(8J!-12J!j' e + QJ e J' e 2 + J' e 3 ) + 3J' c (6J 2 e -12J e J' e + J^ 

- 6J c (4J e 2 + 7J e J' e - J' e 2 ) - 3(4 J 2 + 12 J C J' C - J' c 2 )(2J e + J' e ) 

+ (8J 3 C -12J 2 J' C + 6J C J' 2 + J' C 3 )]. (9) 

If the Hamiltonian (|7|) completely describes the material, we have x HTE {T) 
x(T) or F m = A m for all m, where A m 's are in Table [l]. However the Hamil- 
tonian may only approximately describe the real material. We can estimate 
optimal values for the exchange parameters by minimizing the kth deviation 
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J2( F m-An) 2 . (10) 
\ K m=l 

In the present case, the exchange parameters are J e , J' e , J c and J' c . We 
numerically minimized D 3 with F m 's in eq.(^j) and obtained the optimal 
values as 

J e ~ 480 K, J' e ~ 530 K, 

J c ~ 580 K, ~ 540 K. (11) 

This estimation is the main result of this paper. The result shows that 
inhomogeneity is weak and frustration is strong in this material. Thus the 
observed spin gap originates from frustration. These values yield deviation 
D3 ~ 0.11 and possible reasons of the deviation are discussed in the next 
section. 
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4 Discussion 



We introduced a novel method to estimate exchange parameters. This method 
is fairly general and so applicable to various systems when we compare a 
Hamiltonian to experiment data. In the first step of this method, we con- 
structed an experimental formula for magnetic suscestibility in the power 
series of T /T. To examine the accuracy of this method, we applied it to a 
known Heisenberg model on a simple square lattice. The formula of 

high temperature expansion for this model is already calculated. [03]. Us- 
ing this formula instead of experimental data, we obtained coefficients A TO 's 
in the way of section 2. We confirmed that the coefficients approximately 
repropduce the original formula. 

We applied the method to C&V4O9 to describe it by a two-dimensional 
Heisenberg model (^). The optimal values for exchange parameters are shown 
in eq. (|TT1). The result shows that the strong frustration opens a spin gap 
in this material. The optimal values yield a deviation of D% ~ 0.11. We 
cannot decide now whether or not this deviation is only an error within 
this method itself. We point out other possible origines of the deviation. 
It is possible that the deviation comes from unexpected components, e. g. 
V0 2 , CaV 2 5 , CaV 3 7 , included in samples of CaX^Og. The deviation also 
possibly means that the material includes some degree of freedom which 
cannot be described by a Heisenberg model (0); it might come from itinerant 
effect or effect of degeneracy of atomic orbitals. Anyway what we have done 
is to approximately represent CaV/iOg by the Heisenberg model (^). 

The values of exchange parameters are plausible if similar materials have 
similar values. In CaX^Os we estimated J ~ 600 K by fitting the formula for 
the one-dimensional Heisenberg model to experiment. This value is similar 
to a typical exchange parameter ~ 500 K for CaV 4 9 . In CaV 3 7 , J c ~ J e 
is argued [F20J by using a theoretical result. [[Hj] These results are actually 
similar to our result for CaV 4 Og. 

Using the values of the exchange parameters, we calculated the spin gap 
of the Hamiltonian (|7|) by the numerical diagonalization. The extrapolated 
spin gap A is given as A ~ 170 K. Considering the 10 % deviation of the 
coefficients and errors for extrapolation, this result seems to be consistent 
with the observed value of A ~ 110 K. 
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